Capillary ordering and layering transitions in two-dimensional hard-rod fluids 
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In this article we calculate the surface phase diagram of a two-dimensional hard-rod fluid confined 
between two hard lines. In a first stage we study the semi-infinite system consisting of an isotropic 
fluid in contact with a single hard line. We have found complete wetting by the columnar phase 
at the wall-isotropic fluid interface. When the fluid is confined between two hard walls, capillary 
columnar ordering occurs via a first-order phase transition. For higher chemical potentials the 
system exhibits layering transitions even for very narrow slits (near the one-dimensional limit). 
The theoretical model used was a density-functional theory based on the Fundamental-Measure 
Functional applied to a fluid of hard rectangles in the restricted-orientation approximation (Zwanzig 
model). The results presented here can be checked experimentally in two-dimensional granular media 
made of rods, where vertical motions induced by an external source and excluded volume interactions 
between the grains allow the system to explore those stationary states which entropically maximize 
packing configurations. We claim that some of the surface phenomena found here can be present in 
two-dimensional granular-media fluids. 

PACS numbers: 64.70.Md,61.30.Hn,61.20.Gy 



INTRODUCTION 



The effect of fluid confinement on phase transitions is 
nowadays an active line of scientific research due to the 
direct application of the theoretically predicted surface 
phase diagrams in the nanotechnology industry. Confin- 
ing simple fluids, such as hard [3, Ej or Lennard- Jones 
0| spheres, in a narrow slit geometry, results in a rich 
phase behavior, which has recently been studied in de- 
tail. Liquid crystals confined in nanopores is another 
typical example of systems with important applications 
in the industry of electronic devices. For this reason they 
have been extensively studied using theoretical models 
based on density-functional theory. In particular, capil- 
lary phase transitions exhibited by a nematic fluid con- 
fined between hard walls 0, Hj], or walls favoring a par- 
ticular anchoring [fj], have been predicted. When non- 
uniform liquid-crystal phases, such as the smectic phase, 
are included in the study of confined systems, the result- 
ing surface phase diagrams display a rich phenomenol- 
ogy, which includes wetting transitions, the appearance 
of smectic defects Q, and layering transitions [§]. 

The effect of confinement on two-dimensional fluids is 
also an interesting topic of research. Langmuir monolay- 
ers of lipids on the surface of water have been extensively 
studied in the last hundred years [§], and the discov- 
ery of structures and phase transitions in these systems 
has experienced a dramatic evolution driven by the new 
experimental techniques. Now it is possible to confine 
these two-dimensional systems by external potentials and 
study the influence of the confinement on the molecular 
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packing of surface monolayers. 

Another paradigm of two-dimensional systems where 
the confinement plays an important role is the pack- 
ing structures formed by particles in granular media 
[nj]. The crystallization of a quasi-two-dimensional one- 
component granular-disk fluid has recently been studied 
experimentally [ll[. It was found that the properties 
of the crystal structure obtained (such as packing frac- 
tion, lattice structure, and Lindenman parameter) co- 
incide with their counterparts obtained from MC sim- 
ulations of a hard disk fluid. Recent experiments have 
found non-equilibrium steady states in a vibrated granu- 
lar rod monolayer with tetratic, nematic and smectic cor- 
relations Some of these textures are similar to the 
equilibrium thermodynamic states of two-dimensional 
anisotropic fluids resulting from density-functional cal- 
culations and MC simulations [T3|. It was shown by 
several authors that the inherent states of some frozen 
granular systems can be described by equilibrium sta- 
tistical mechanics [15]. Also, an experimental test of 
the thermodynamic approach to granular media has re- 
cently been carried out [16]. Confining two-dimensional 
granular rods in different geometries (circular, rectangu- 
lar, etc.) results in spontaneous formation of patterns, 
with different orientationally ordered textures and de- 
fects next to the container [TtJ • The authors of Ref. [l8| 
have carried out MC simulations of a confined hard disk 
fluid. They have found that the crystal phase fails to nu- 
cleate due to formation of smectic bands when the system 
is confined [18]. It would be interesting to device an ex- 
periment with confined granular disks with the aim of 
comparing the properties of the non-uniform stationary 
states with those obtained from the statistical mechanics 
applied to the hard disk fluid. 

The main purpose of this article is the study of a con- 
fined two-dimensional hard-rod fluid. We are interested 
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in the calculation of the surface phase diagram of a hard 
rectangle (HR) fluid confined by a single or two hard 
lines. We can think on a HR fluid as an experimental re- 
alization of a system of hard cylinders confined between 
two plates at a distance less than twice the cylinder di- 
ameter. We suggest that some of the surface phase tran- 
sitions obtained here by applying the density functional 
formalism to a confined two-dimensional HR fluid should 
be similar to the steady states of confined granular rods. 
Some experiments are required to verify this hypothesis. 

The paper is organized as follows. In Sec. II we 
present the theoretical model: the fundamental-measure 
density functional applied to a HR fluid in the restricted- 
orientation approximation. This section is divided into 
two subsections. In the first the model is particularized 
to the study of the bulk phases, while in the second part 
the theoretical expressions used in the calculations of the 
thermodynamic and structural properties of the inter- 
faces are presented. The results are presented in Sec. 
III. First we study the bulk phase diagram of a HR fluid 
with aspect ratio equal to 3, and then the resulting sur- 
face phase diagrams of a single wall-HR fluid interface 
and of the fluid confined between two hard lines are pre- 
sented. Some conclusions are drawn in Sec. IV. 
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and where the symbol * stands for convolution, i.e., 
p v * ui*' = J v dr' p t ,{r')uj < f'\v — r'). The weights 
are the characteristic functions whose volume integrals 
constitute the fundamental measures of a single particle 
(the edge lengths and surface area). They are defined as 
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where = a + (L — a)6fj, v , with L and a the length and 
width of the rectangle and the Kronecker function, 
while S(x) and Q(x) are the Dirac delta and Heaviside 
functions, respectively. 



II. THEORETICAL MODEL 

In this section we introduce the theoretical model used 
in the calculations of the bulk and interface equilibrium 
phases. To study highly inhomogeneous phases such as 
those resulting from the confinement of a fluid in a nar- 
row slit geometry or the solid phase with a high packing 
fraction, we have used the Fundamental-Measure Theory 
(FMT) applied to an anisotropic fluid of hard rectangles. 
It is well known that this formalism presents a great ad- 
vantage over other techniques when dealing with highly 
inhomogeneous phases, and that this is mainly due to the 
fact that a basic requirement to construct the FMT den- 
sity functional is that it conform with the dimensional 
cross-over criterium [13, [23]. To implement the calcu- 
lations we have used the restricted-orientation approxi- 
mation, where the axes of the rectangles are restricted 
to align only along the coordinate axes x or y. Thus, 
the whole system is described in terms of density profiles 
p v (r) (v = x,y). 

While the ideal part of the free energy density in re- 
duced thermal units has the exact form 

Mr) = EM')M')-i]> (i) 

V 

the FMT interaction part of the 2D HR fluid is approxi- 
mated [23] by 

$ cxc (r) - -no(r) In [1 - n 2 (r)] + n ^ n ^\ (2) 



A. The bulk phases 

To calculate the bulk phase diagram we need to mini- 
mize the Helmholtz free energy functional /3J r [{p l , (r)}] = 
J dr[$id(r) + <£> oxc (r)] with respect to the density pro- 
files p v (r). These density profiles have the symmetries 
corresponding to the equilibrium phases, which can be 
the isotropic or nematic fluids, the smectic phase (with 
particles arranged in layers with their long axes pointing 
perpendicular to the layers), the columnar phase (with 
long axes parallel to the layers), plastic solid (particles 
located at the nodes of the square grid with averaged ori- 
entational order parameter over the cell equal to zero), 
and oriented solid (with both translational and orienta- 
tional order) . To take proper account of all these possible 
symmetries, we have used a Fourier-series expansion of 
the density profiles: 

N 

Pu{v) = pQX v ^2 a( kuk 2 COS (H X ) COS (<l2y)' ( 8 ) 
k=(0,0) 

where we defined k = (fci,fc 2 ) [with N = iV(l,l)], 
qi = 2wki/d x , and 172 = 2wk2/d y are the wave vector 
components parallel to x and y axes respectively, and d x , 
d y are the periods of the rectangular grid along these di- 
rections. a£ k are the Fourier amplitudes of the density 

profile of the species v with the constraint olq \ = 1. po is 
the average of the local density over the cell, while x v is 
the cell-averaged occupancy probability of species v. The 
Fourier series is truncated at that value N which guar- 
antees that ocjlpjf < 10~ 7 . With this parametrization the 
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weighted density can be calculated explicity as 

n a (r) = p ^x v a^ k2 ^i a) 0i-)cos(qix)cos(q2y), (9) 

where 6jl a \h) are the Fourier transforms of the corre- 
sponding weights, which have the form 

4°)(k) = xo(«io5/2) X o(ga«7»/2) ) (10) 
u;i lx \k) = <xi(<?i^/2)xoM/2), (11) 

4^)(k) = a*xo(«io5/2)xi(gi l o*/2), (12) 
^ 2 )(k) = a X i(gi^/2) Xl (g 2 ^/2), (13) 

Here a = Lcr is the surface area of the particle, and 
Xo{ x ) = cos a;, Xi( x ) — sm(x)/x. We have selected the 
orientational director parallel to y. Thus, the equilib- 
rium smectic (columnar) phase should be found by min- 
imizing the free energy with respect to the Fourier am- 
plitudes a^l (cKfco): the smectic (columnar) period d y 
(d x ) and the order parameter Qn £ [— 1,1] [related to 
the x u 's through the relations x\\ t ± = (1 ± Qn)/2 where 
the symbols ||,_L stand for particle alignment along y 
and x respectively]. For uniform phases [a^\ = 
V(fci,&2) 7^ (0,0)] Qn coincides with the nematic order 
parameter. The solid phase is to be found by minimiz- 
ing the free energy with respect to all the Fourier ampli- 
tudes ot.fa k , the crystal periods d x and d y , and the order 
parameter Qn in the case of an orientationally ordered 
solid. To measure the packing structure and the orienta- 
tional order of the bulk phases we use the local density 
and the order parameter profiles, p(r) — J2vP"( r )> anc ^ 
Q(r) = [p y (r) - p x {r)} / p{r) respectively. 



favored homeotropic alignment usually considered in sev- 
eral studies of three-dimensional liquid crystals confined 
by a single or two walls (in particular that of Ref. [H]). 

The one-dimensional equilibrium density profiles p v (x) 
were found by minimizing the excess surface free energy 
per unit length 

7 =y dT |£M + p_^^ (x) [^_ K(x)] | ; ( i6) 

where = {k B T)- 1 , ®(x) = $ id (x) + <f> cxc (x), and p v 
are the chemical potentials of species v fixed at the bulk 
fluid-phase value at infinite distance from the wall, while 
P is the fluid pressure. The chemical potential of the 
bulk fluid phase is calculated, as usual, as p = Yl v 
with x v the molar fractions of species v. If the bulk phase 
is an isotropic fluid then x v — 1/2, and p v = p, W. 

To measure the degree of interfacial order, we will 
use the adsorption of the density profile, defined as 
r = J dx [p{x) — p(oo)], and the order parameter profile 
Q(x). 

The expression l|16p coincides with the definition of the 
surface tension of the wall-fluid interface for the semi- 
infinite case, which is approximately equal to half the 
excess surface free-energy for the slit geometry when the 
wall distance H is large enough to accommodate both 
interfaces. 

To minimize the functional given by (fT6|) , we have dis- 
cretized space in the x direction and minimize 7 with 
respect to p v (xi) (xi 6 [xq^xn]) using the conjugate- 
gradient algorithm. 

III. RESULTS 



B. The interfacial phases 

As we want to study the hard wall-fluid interface or the 
HR fluid confined in a slit geometry, we have introduced 
the following external potential: 

( 00, x < a x v /2 
Vu{x) = { (14) 
[0, x > o£/2, 

for the semi-infinite system, and 

!oo, x<a^/2 and x>H — af,/2 
(15) 
0, al/2<x<H-cjll2, 

for the slit geometry, where H is the slit width, and 
the normal to the wall was selected in the x direction. 
Note that this external potential represents a hard line 
which excludes the center of mass of particles at distances 
less than their contact distances with the wall. In this 
sense we can say that the external potential favors par- 
allel alignment at the wall. This is in contrast with the 



In this section we present the main results obtained 
from the application of the theoretical model just de- 
scribed to the study of surface properties of a 2D HR 
fluid. Particles were chosen to have aspect ratio k 
I. n = 3. This aspect ratio is chosen because one of the 
aims of the present work is the study of layered phases 
confined by one or two walls. As we will show bellow for 
k — 3 the stable phase is the columnar layered phase. 

In the first subsection we will summarize the results 
obtained in the calculation of the bulk phase diagram of 
this system, while in the second subsection we will focus 
on the study of the surface phase diagram. 

A. Bulk phase diagram 

We have minimized the free energy density of the HR 
fluid, defined as $ = V^ 1 J v dr [$ id (r) + $ oxc (r)], with 
respect to the Fourier amplitudes, periods, and mean oc- 
cupancy probability, as described in detail in Sec. Ill Al 
The results are plotted in Fig. [U where the free-energy 
densities of all the stable and metastable phases found 
are plotted as a function of the packing fraction r\ = poa. 



We have found, apart from the usual isotropic (I) and ne- 
matic (N) phases, two different smectic phases (Smi, and 
Sma), a plastic solid (PS), perfectly oriented solid (OS), 
and finally the columnar phase (C), which is the stable 
one in the whole range of packing fractions explored. 
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Figure 1: The rescaled free-energy density $* = $ + 2.9875 — 
5.850177 is plotted against the mean packing fraction for all 
the stable and metastable phases found. These are: isotropic 
(dashed line), nematic (dotted line), smectic- 1 and smectic- 
2 (dotted and dashed lines), plastic solid (dashed line la- 
belled as PS), while the perfectly oriented solid and the 
columnar phases (labelled in the figure as OS and C respec- 
tively) are plotted with solid lines. The open circle indi- 
cates the isotropic-nematic bifurcation point; the open square, 
the isotropic-plastic solid bifurcation point; and the solid cir- 
cles represent the coexisting packing fractions at isotropic- 
columnar phase coexistence. 

The coupling between the spatial and orientational de- 
grees of freedom of the particles results in the presence of 
phases (stable or metastable) with different symmetries. 
In Fig. ([2]) we have sketched some of the particle con- 
figurations corresponding to phases with columnar (a), 
smectic- 1 (b), smectic-2 (c), and plastic solid (d) sym- 
metries found from the numerical minimization of the 
density functional. The directions of spatial periodicities 
of each phase have been depicted in the figure. 

In Fig. 02 (a) we have plotted the density and order- 
parameter profiles of the coexisting columnar phase. The 
columnar phase is orientationally ordered in the y direc- 
tion with the long rectangle axis pointing along the y axis, 
while the periodicity of both density and order parame- 
ter profiles (which are in phase) is along the x direction 
[see Fig. [3] (a)]. The mean coexistence packing fractions 
of the I and C phases are rji = 0.57058 and rye = 0.60310, 
respectively while the period of the C phase, in units of 
the HR width, was found to be d x /a = 1.20102. In Fig. 
[3] (b) we have plotted the order parameter Qn, and the 
period of the columnar phase as a function of the mean 
packing fraction. 
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Figure 2: Sketch of particle configurations corresponding to 
different phases: columnar (a), smectic-1 (b), smectic-2 (c), 
and plastic solid (d) phases. The direction of spatial period- 
icities are labeled in the figure. 



To compare the different packings of HR particles in 
the metastable phases (found as the local minima of the 
free energy density) for a fixed mean packing fraction 
i] = 0.7, we have plotted the density and order-parameter 
profiles of the Srni >2 [Fig. g](a) and (b)], and PS and OS 
[Fig. [5] (a)-(c)] phases. As can be seen from Fig. [4] 
(a), the density profile of the Smi phase has two max- 
ima per period. The less pronounced maxima, located 
at the interstitials, reflect the high population of parti- 
cles with long axes oriented parallel to the smectic layers 
[see the sketched particle configurations in Fig. [2] (b)]. 
This alignment is also shown in the order-parameter pro- 
file, which reaches high negative values at the intersti- 
tial positions. This phase bears a strong resemblance 
to the findings of Refs. [2l[ and [13] where the parti- 
cle equilibrium configurations in the 3D smectic phases 
show the same pattern. As a consequence of this (al- 
ternating population of particles aligned perpendicular 
-sharpest peak in the density profile- and parallel to the 
layers) , the smectic period in units of the particle length 
is dy/L — 1.53025, higher than the smectic period of the 
Sm2 phase (dy/L = 1.17935). The density and order- 
parameter profiles of the Sm2 are shown in Fig. [4] (b). 
As can be seen from the figure, these profiles reflect the 
usual packing in smectics, characterized by a single den- 
sity peak with vanishingly small population of particles 
in the interstitials, while the order parameter reaches its 
maximum value at the position of the smectic layers [see 
Fig. [2] (c) for the sketched particle configurations] . 

The density and order parameter profiles of the PS 
phase with mean packing fraction equal to 0.7 are plot- 
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Figure 3: (a): density p{x) (solid line) and order parameter 
Q(x) (dashed line) profiles of the columnar phase at coex- 
istence with the isotropic phase, (b): order parameter Qn 
and period of the columnar phase against the mean packing 
fraction. 




Figure 4: Density (solid line) and order parameter (dashed 
line) profiles of the smectic-1 (a) and smectic-2 (b) metastable 
phases for a value of mean packing fraction fixed at 0.7. 



ted in Fig. [5] (a) and (b). The plastic solid has the same 
periodicity in the x and y direction, i.e. d x = d y = d, 
and the order parameter averaged over the unit cell is 
strictly equal to zero. As we can see from Fig. [5] (b), 
while the order parameter at the nodes of the square 
lattice is equal to zero, it reaches positive (negative) val- 
ues at the (±0.5,0) [(0, ±0.5)] positions along the sides 
of the cell (the same solution with the x and y direc- 
tions interchanged was found in the minimization of the 
free energy) . Finally, the density profile of the perfectly 
aligned two-dimensional solid is plotted in Fig. [5] (c). 

Although the phases described above are metastable 
with respect to the columnar phase, they can be stabi- 
lized for different values of the particle aspect ratio. A 
detailed study of the complete phase diagram, necessary 
to elucidate this point, is a work in progress. 

We now proceed to make a comparison between the 
results for the 2D Zwanzig model with k = 3 obtained 



above and those for hard parallelepipeds with restricted 
orientations and the same value of k [23]. This compari- 
son will show the differences in phase behaviour between 
three and two dimensions as predicted by Fundamental- 
Measure Theory (which, as already pointed out, con- 
forms with the dimensional crossover criterion). As 
shown in Ref. [23| , hard parallelepipeds exhibit a second- 
order phase transition between isotropic and plastic solid 
phases. As density increases the system goes to a discotic 
smectic phase (confirmed by simulations) via a first-order 
phase transition, which in turn discontinuously changes 
to a columnar phase and then to an oriented solid. By 
contrast, the present model shows that, in two dimen- 
sions, the isotropic phase exhibits a first-order transition 
to a columnar phase that is stable until very high pack- 
ing fractions (more stable that plastic, oriented solid and 
different smectic phases). As a consequence, one expects 
that the corresponding surface phase diagrams be also 
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Figure 5: Density (a) and order-parameter (b) profiles of the 
plastic solid phase, (c): Density profile of the perfectly ori- 
ented solid. 



different. 



B. Surface phase diagram 

In this section we deal with surface phenomena. In the 
first part we will concentrate on the semi-infinite wall- 
isotropic interface of a HR fluid, while in the second part 
we will focus on the slit geometry. We will demostrate 
the presence of complete wetting, capillary ordering and 
layering transitions in the confined two dimensional hard 
rod fluid. For a detailed discussion on general grounds of 
the phase behavior and critical phenomena of a confined 
by a single wall fluid see Ref. [24(. 



1. The wall- fluid interface. 

The interaction between the isotropic fluid phase and a 
hard wall was studied by calculating the one-dimensional 
equilibrium density p{x) and order-parameter Q(x) pro- 
files through the excess surface free-energy minimization 
[see Eq. (|16|l]. The chemical potential fi of the fluid phase 
at infinite distance from the wall was varied within the 
range of isotropic-phase stability, i.e. [i G [— oo,/i ] (/u 
being the value at which the I-C phase transition occurs). 
It is well known that the presence of a hard wall in a sys- 
tem of elongated particles induces parallel alignment of 
the particle axes with respect to the wall jlEIlfij]- This 
preferential alignment is a result of the entropic depletion 
effect. In the parallel configuration, the centers of mass 
of the particles are much closer to the wall, so the gain 
in volume per particle is larger and, as a consequence, 
the configurational entropy of the system is maximized. 
This effect is responsible for the occurrence of a biaxial 
nematic phase which breaks the orientational symmetry 
in a three-dimensional nematic fluid The same deple- 
tion mechanism is at work in 2D, as we will show below. 

The results from the minimization are shown in Figs. [6] 
(a) and (b) for an undersaturation of /3A/i = —1.1 x 10~ 4 . 
As we can see from the figure, the density and order- 
parameter profiles indicate columnar order near the wall, 
which propagates several columnar periods into the fluid 
phase. The wall-fluid interaction enhances the orienta- 
tional order near the surface and the adsorption of par- 
ticles, creating a structured layer with columnar-phase 
symmetry which grows in width with increasing chem- 
ical potential and diverges at fi = /xq- Thus, complete 
wetting by a columnar phase occurs at the wall-isotropic 
interface. This result is shown in Fig. [3(a) where the ex- 
cess surface free-energy 7 and the adsorption coefficient T 
are plotted against A/j, = /i — [j,q. As we can see, T grows 
continuously, ultimately diverging logarithmically with 
A/z (see inset of figure) . The excess surface free energy 7 
has a maximum, and at this point the adsorption passes 
through zero. This result is directly related to the interfa- 
cial Gibbs-Duhem equation, T = —d^/dp, which relates 
the adsorption coefficient with the first derivative of the 
excess surface free-energy with respect to bulk chemical 
potential. At the excess surface energy is equal to the 
wall-isotropic surface tension 7wi, which is in turn equal 
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to the sum of wall-columnar and columnar-isotropic sur- 
face tensions, 7(^0) = 7wi = 7wc + 7ci (the Young's 
equation for complete wetting). 
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Figure 6: Density (a) and order parameter (b) profiles of the 
wall-isotropic fluid interface. The undersaturation is fixed to 
fiAfj, = -1.1 x 10~ 4 . 

We have carried out a logarithmic fit of the adsorp- 
tion coefficient with respect to undersaturation PA/i = 
(3 ({j, — (Xo), and we find that Ta = Ti +r 2 In [/3|A/x|], with 
ri = 0.02387 and t 2 = -0.03396. Then, integrating the 
interfacial Gibbs-Duhem relation T — —dj/dfj., we find 
the expression 



(3jcr « ^7wi<t 



n 



r 2 (ln(/3|A At |)-l)]/3A M , (17) 



which approximates the excess surface free energy near 
complete wetting. The above expression is plotted 
against /3A/i in Fig. [7] (b) , where the results from direct 
calculation of P^ya, using the equilibrium density profiles 
obtained, are also plotted. As we can see the agreement 
is excellent even for relatively high values of undersatu- 
ration. 
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Figure 7: (a): excess surface free energy (solid line) and ad- 
sorption coefficient (dashed line), in reduced units, against 
/3A^i. The inset shows Tcr vs. /3A/i in logarithmic scale, (b): 
excess surface free energy vs. f3fi in the neighborhood of zero 
undersaturation. The open circles show the values obtained 
from the numerical minimization, while the solid line repre- 
sents the analytic curve obtained by integrating the interfacial 
Gibbs-Duhem relation with the fitted logarithmic dependence 
of the adsorption coefficient (see text). The solid circle shows 
the value of the W-I surface tension /37wic = 0.13822 



To calculate the structural and thermodynamic prop- 
erties of the columnar-isotropic interface, we have imple- 
mented a numerical scheme already used in Ref. |27l |. 
consisting of minimizing the surface excess free energy 7 
in a box of width h containing a stripe of a few colum- 
nar layers surrounded by isotropic material with periodic 
boundary conditions, h is chosen such that the density 
profiles can easily accommodate the two interfaces and 
go to the coexistence fluid density at the periodic bound- 
ary. A typical result from this calculation is plotted in 
Fig. [8] (a) and (b) for the density and order-parameter 
profiles, respectively. Thus, the I-C interfacial tension 
can be calculated as half the excess surface free energy 
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resulting from the minimization. We have found a value 
of /3 7ic cr = 0.00672. 
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Figure 8: Density (a) and order-parameter (b) profiles of a 
numerical box containing two isotropic-columnar interfaces. 

Finally, to verify that Young's law for complete wet- 
ting holds, we need to calculate the surface tension of 
the wall-columnar interface. To construct density pro- 
files compatible with this semi-infinite interface, one has 
to establish a boundary, at the side of the computational 
box opposite to the wall, and place, beyond the bound- 
ary and into the bulk, a periodically structured profile, 
choosing the phase (i.e. the value of the profile at the 
boundary) arbitrarily within the bulk period. Although 
this recipe can in principle be implemented, we have cho- 
sen to fix bulk I-C coexistence conditions in a confined 
columnar phase and calculate the density profile of the 
system confined between two walls; the separation be- 
tween the walls was chosen large enough so that the ef- 
fects of having a finite interface penetration length caused 
by the presence of the confined external potential can be 
neglected. Also, in order to ensure that commensurabil- 
ity effects can be ignored, the distance between the walls 



was set to a (large) integer number of equilibrium periods 
of the columnar phase. The results from these calcula- 
tions are plotted in Fig. [9] (a) and (b). The W-C surface 
tension calculated as half the value of the excess surface 
free energy results in /37wcc = 0.13150, compatible with 
Young's law in conditions of complete wetting of the W-I 
interface by the columnar phase. 
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Figure 9: The density (a) and order parameter (b) profiles 
of two wall-columnar interfaces 



2. 



This section is devoted to a study of the effect of con- 
finement of a 2D HR fluid on the thermodynamic and 
structural properties of the fluid. In particular, we are 
interested in the enhancement of the orientational and 
layering ordering due to confinement, and the commen- 
surability effects exhibited by a layered phase sandwiched 
between two hard walls at a distance that may or may 
not commensurate with the period of the bulk columnar 
phase. It is well known that, under certain circumstances 
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(related to the nature of the fluid-fluid and surface-fluid 
interactions), a fluid inside a pore can exhibit capillary 
first-order phase transitions between two different phases 
at a chemical potential below the bulk coexistence value. 
An example of this phenomenon is the recently studied 
capillary nematization [3] and smectization || of a liquid 
crystal fluid inside a pore. The bulk condensed phase 
may have uniform (nematic) or nonuniform density pro- 
files. For the latter case, capillary layering transitions be- 
tween interfacial phases with different number of smectic 
layers [|| can also be found. Here we will show that these 
capillary and layering phase transitions are not unique to 
3D system. They are also present in 2D anisotropic fluids 
which can stabilize layered phases with different spatial 
symmetries, such as the columnar phase. 

With a view to finding the effects of confinement on 
columnar ordering in a HR fluid, we have minimized the 
excess surface free energy with respect to the density pro- 
file for the particular case of HR's with k = 3. The fluid 
is confined by two hard walls at a distance H/ a = 30 (in 
units of the particle width). As already pointed out, hard 
walls favor alignment parallel to the wall, as well as ad- 
sorption of particles at both surfaces (density and order 
parameters at contact are much higher than their bulk 
values). This coupled translational-orientational order- 
ing near the surfaces propagates into the fluid, creating 
columnar ordering. We have found that for low values 
of the chemical potential of the bath the density profile 
is structureless (except just at the wall contact), similar 
to the bulk isotropic phase. Increasing the chemical po- 
tential several damped columnar peaks appear near the 
wall in a continuous fashion, i.e. with their heights in- 
creasing continuously. At some value of the chemical po- 
tential, the system exhibits a first-order phase transition 
between a phase with highly damped columnar peaks to 
a new phase with much stronger columnar ordering even 
at the center of the pore. The typical density and order- 
parameter profiles of both interfacial phases are shown in 
Fig. [I0](a)-(d). Although the less-ordered phase exhibits 
strong oscillations in both density and order-parameter 
profiles, the peak amplitudes are damped into the pore 
faster than those of the higher ordered phase. We will 
take the convention to call the first 'isotropic', and the 
second 'columnar' surface phases. This convention is jus- 
tified by the fact that, just before the transition described 
above, columnar ordering increases continuously, start- 
ing from an isotropic-like density profile, as the chemi- 
cal potential is increased. Thus we cannot trace out a 
definite boundary (a value for /i below that correspond- 
ing to first order phase transition) below or above which 
the profile inside the pore can be considered isotropic or 
columnar. Only the first order phase transition described 
above can really distinguish two different surface phases, 
one of them less ordered (following our convention, the 
isotropic phase) than the other (the columnar phase). As 
we can see in the figure, the latter has 25 columnar peaks. 

The transition point is calculated from the discontinu- 
ity in the first derivative of the excess surface free en- 
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Figure 10: Isotropic (a)-(b) and columnar (c)-(d) phases that 
coexist at the same chemical potential bellow no. (a), (c): 
Density profiles, (b), (d): Order-parameter profiles. 



ergy with respect to the bulk packing fraction i]. The 
corresponding plot is shown in Fig Qj] (a) . At this point 
the adsorption coefficient jumps discontinuously from the 
less- (the damped columnar) to the higher-ordered phase 
[see Fig. HH(b)] . 

This surface transition point is located below the bulk 
isotropic-columnar phase transition [see Fig. Qj] (a)], 
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Figure 11: Excess surface free energy (a) and adsorption coef- 
ficient (b) against packing fraction of the bulk isotropic fluid. 
In the figure at top, the solid circle represents the transition 
point between both interfacial phases, while open square in- 
dicates the point corresponding to the bulk coexistence value 
for isotropic and columnar phases. 



showing the presence of columnar-order enhancement in 
the pore. On further increasing the chemical potential 
up to a sufficiently high value (above fio), we find a first- 
order layering transition between two interfacial colum- 
nar phases which differ by just a single columnar layer. 
The behavior of the excess surface free energy and the 
adsorption coefficient is similar to that shown in Fig. Ql] 
(a) and (b) . Alternatively we can find the transition from 
n-lton columnar layers by fixing the chemical potential 
and increasing the pore width H. 

The two surface phase transitions described above, 
namely first-order capillary I-C ordering and (n — 1)~ 
n layering transition, are connected in the fi — H surface 
phase diagram through the peculiar structure shown in 
Fig.[H 

The parabola below the bulk transition line corre- 
sponds to first-order transition lines separating regions 
of stability of the isotropic and the columnar interfacial 
phases, while the straight lines indicate layering transi- 
tions. Increasing the chemical potential from low values 
to those corresponding to the parabola, the density pro- 
files always change continously from a structureless to 
damped columnar density profile. Both types of transi- 
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Figure 12: ft — H surface phase diagram showing first-order 
capillary columnar ordering and layering transitions. The 
pore width covers a range which goes from 23 to 25 colum- 
nar layers, as labeled in the figure. The open circles indicate 
the cases chosen for calculations, while the solid line is a cu- 
bic spline interpolation. The horizontal dashed line shows the 
value of the bulk chemical potential at the isotropic-columnar 
phase coexistence. 



tions (the isotropic-columnar and n — 1-n layering transi- 
tions) coalesce in triple points, two of which are shown in 
Fig. [121 At the triple points an isotropic and two colum- 
nar interfacial phases with n — 1 and n layers coexist in 
equilibrium. The set of connected of Fig. [12] are similar 
to those found in MC simulations of the confined hard- 
sphere fluid [3] • In this work the authors have shown the 
existence of capillary freezing of the HS fluid, confined 
in the slit geometry, for chemical potential values below 
the bulk freezing transition. The transitions lines in the 
fi— H surface phase diagram follow the same topology of 
connected set of parabolas as found in our system. 

Some of the topological features of this surface phase 
diagram can be elucidated from the Clausius-Clapeyron 
equation as applied to the interfacial coexistence lines. 
The excess surface free energy ^(fi, H) along coexistence 
is a function of two variables, the chemical potential fi, 
and the pore width H. Thus, infinitesimal changes in 
these variables along the coexistence curve are related 
through the equation 



£?7 

dfi 



dfi 



H 



~dH 



dH = 0, (18) 



where the coexisting condition j a = 7/3 (for 
a, f3 =I,C n _i,C n ) was used, and Au = u a — up for any 
function u. Using the interfacial Gibbs-Duhem equation 
dj/dfi = — r and the definition of the solvation force 
/ = —dj/dH, we arrive at 



dfi 
dH 



AT' 



(19) 
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which relates the first derivative of the chemical potential 
with respect to the pore width with changes in the sol- 
vation force and in the adsorption coefficient at the tran- 
sition point. The negative slope of the layering curves is 
a direct result of Eq. lfT9|) . as the increment in the ad- 
sorption is always positive for the (n — 1) — > n layering 
transition, while the change in the solvation force is also 
positive (the latter can be interpreted as an increment 
with respect to the bulk of the excess surface pressure, 
which is obviously larger for the phase with n layers). 
For values of the pore width that commensurate with an 
integer number of columnar periods of the bulk columnar 
phase, the solvation force becomes zero and we get a min- 
imum in the I-C capillary transition curve (see Fig. 112)1 . 
At each side of the minimum the solvation force change 
the sign to positive (left side) or negative (right side) as 
we compress or expand the film, respectively, while the 
change in adsorption remains positive. 

The Kelvin equation for capillary condensation in a slit 
geometry relates the undersaturation in chemical poten- 
tial with pore width H as 



A/i = p,(H) — hq = — 



(Pa - Pp)H' 



(20) 



where p a and p$ are the bulk coexisting densities of 
phases a and (3 (a being the condensed phase), while 
7 Q( 3 is the surface tension of the corresponding interface. 
It was assumed that complete wetting by the a phase 
occurs at the W-/3 interface. For a detailed discussion of 
the Kelvin equation in the context of liquid crystal phase 
transitions see Ref. (28|. Applying this equation using 
H/a = 28.88 (the location of the minimum in the \i — H 
phase diagram of Fig. [12)) , we obtain an undersaturation 
@Afi = -0.0429, while its real value is (3Afi = -0.1255. 
In the derivation of the Kelvin equation, deviations from 
the bulk structure of the density profile arising from 
the confinement by the external potential are neglected. 
Also, the elastic energy resulting from the compression or 
expansion of a layered phase confined between two walls 
is not taken into account. These effects might be respon- 
sible for the differences found between our calculations 
and the estimation based on the Kelvin equation. We 
have checked that the sequence of minima in the [i — H 
phase diagram tends to [io as H — > oo, a result predicted 
by Eq. (|20l) . 

Refs. [ij] and @ showed that the capillary nemati- 
zation line of the confined liquid-crystal fluid ends in a 
critical point for small values of the pore width. In or- 
der to study how the topology of the surface phase dia- 
gram changes in the regime of small pore widths, we have 
carried out the corresponding calculations of interfacial 
structure. We have found that the I-C capillary ordering 
transition changes at some particular value of H (near 
its maximum undersaturation represented by the mini- 
mum in the I-C interface coexisting curve) from first to 
second order. For lower values of H two critical points 
emerge from this single point, the distance between them 
increasing. In Fig. [13] one of these scenarios is shown. 



As we can see, there is a range of values of H (near 
the triple points) where the first-order capillary order- 
ing transitions are still present but, between the critical 
points belonging to different layering branches, columnar 
ordering grows continously from the isotropic (damped 
columnar interfacial phase) to a highly-ordered colum- 
nar phase. Layering transitions are always present even 
for very small H, as will be shown below. An interest- 
ing feature of this phase diagram is that the location of 
the triple points moves above the bulk coexistence value 
/io- This indicates that the interfacial columnar phase 
just below the triple points can be unstable for values of 
chemical potentials corresponding to those of columnar- 
phase stability at bulk (similar to the capillary evapora- 
tion of the confined fluid) . For wide enough slits (those 
for which the parabolas are connected) the triple points 
are practically located at as can be observed from 
Fig. m 
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Figure 13: \i—H surface phase diagram showing critical points 
(rilled circles). Number of columnar layers are indicated as 
subscripts. 

For even smaller values of H, only layering transitions 
remain; these end in critical points located above /io, as 
Fig. [TD shows. 

When the width H is such that the pore can only ac- 
commodate one particle with its long axis perpendicu- 
lar to the wall (or not more than four or three particles 
aligned parallel to the wall) the system is near the one- 
dimensional limit. It is known that hard-core systems in 
this limit do not exhibit first-order phase transitions, but 
even for very narrow slits we can still find first-order tran- 
sitions at which the density profile experiences an abrupt 
change inside the pore. In Fig. [15] (a) and (b) we show 
two coexisting density profiles corresponding to oversat- 
urations, /3A/z = 0.51760 and /3A/i = 0.72836, and pore 
widths H/a = 4.32 and H/a = 3.14 in (a) and (b), re- 
spectively. The fluid inside the pore undergoes a phase 
transition, which dramatically changes the structure of 
the interfacial density profiles by increasing the heigh of 
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four [Fig. flBTa)] or three [(b)] density peaks inside the 
pore. 



IV. CONCLUSIONS 

In this article we have shown that 2D fluids composed 
of anisotropic particles interacting via hard-core repul- 
sion and confined in a slit geometry exhibit a complex 
and rich interfacial phase behavior. Apart from first- 
order capillary columnar ordering, we have also found 
layering transitions in this system. These results are sim- 
ilar to those found in 3D liquid-crystal fluids confined in 
a pore, where capillary smectization and layering phe- 
nomena were also found [8j. In view of these similarities, 
we can extract the conclusion that, independent of the 
system dimensionality and the peculiarities of the layered 
phases, either smectic or columnar, if the fluid- wall inter- 
action enhances layered interface ordering (homeotropic 
in case of smectic phases, and the entropically favored 
parallel alignment for the columnar phase), compatible 
with the equilibrium bulk phase, then the confined fluid 
exhibits the interfacial phase transitions described above. 

In this study we have used as a model a hard-rectangle 
fluid, and the density and the order-parameter pro- 
files were calculated by minimizing the excess surface 
free-energy functional resulting from the Fundamental- 
Measure Theory applied to the two-dimensional Zwanzig 
model. The orientational degrees of freedom were 
discretized, in order to take advantage of having a 
free-energy functional which reduces to the exact one- 
dimensional functional when the density profile is con- 
strained to lie along a line. This property is crucial 
to study strongly confined fluids (as is the case in this 
study), in particular when the pore width has only a few 
particle diameters in width. 

As already pointed out in Sec. [J some experiments 
had shown profound similarities between particle con- 
figurations obtained as stationary states of systems of 
anisotropic grains and those corresponding to the equi- 
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Figure 15: (a): density profiles of two coexisting phases 
(shown with solid and dashed lines) at /3Afi = 0.5176. The 
pore width is H/o = 4.32. (b): same as in (a) but for a pore 
with H/o = 3.14 and for /3A/J, = 0.72836. 



librium states obtained by density functional minimiza- 
tion [l3| ■ These similarities can be explained by applying 
a maximum-entropy principle on granular collections of 
particles, i.e. for a fixed packing fraction, externally- 
induced vibrational motion forces the system to explore 
those stationary states which maximize the configura- 
tional entropy (since the grains cannot overlap). Of 
course, equilibrium statistical mechanics is unable to pro- 
pose an equation of state for granular matter, but it could 
be possible to predict that granular matter composed of 
anisotropic particles and confined between parallel walls 
may support a stationary texture consisting of layers of 
particles oriented parallel to the wall. The manner in 
which the grain orientations propagate into the container 
would depend on the average packing fraction and on 
the frequency of the external force. Only at this qualita- 
tive level can we give some insight into possible complete 
wetting phenomena and capillary ordering in granular 
rod fluids confined between two horizontal plates at a 
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distance slightly larger than the particle dimensions in 
the vertical direction (thus simulating a two-dimensional 
system), and also confined by one or two vertical planes 
(these playing the role of hard walls) . 

Some calculations (not shown here) on the 2D HR fluid 
show that, for different aspect ratios, 2D smectic and 
crystal phases can be stable over some range of pack- 
ing fractions. It would be interesting to explore whether 
confinement suppresses or enhances bulk ordering, and 
to study the changes in the surface phase diagram when 
phases of different symmetries are included. Work along 
this direction is currently in progress. 
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